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Abstract 



>• I We introduce a new family of nonperiodic tilings, based on a substitution rule 



that generalizes the pinwheel tiling of Conway and Radin. In each tiling the tiles 

^ ■ are similar to a single triangular prototile. In a countable number of cases, the 

^ . tiles appear in a finite number of sizes and an infinite number of orientations. 

f^^ ! These tilings generally do not meet full-edge to full-edge, but can be forced through 

On I local matching rules. In a countable number of cases, the tiles appear in a finite 

^ ■ number of orientations but an infinite number of sizes, all within a set range, while 

c^ ■ in an uncountable number of cases both the number of sizes and the number of 

p . orientations is infinite. 
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§1. Introduction 

We introduce a new family of nonperiodic tilings, indexed by a continuous 
parameter. In each tiling the tiles are similar to a single triangular prototile. In a 
countable number of cases, the tiles appear in a finite number of sizes and an infinite 
number of orientations. In a countable number of cases, the tiles appear in a finite 
number of orientations but an infinite number of sizes, all within a set range. In 
one case both the number of sizes and orientations is finite, while in an uncountable 
number of cases both the number of sizes and the number of orientations is infinite. 
A piece of a tiling with two sizes and an infinite number of orientations is shown in 
Figure 1. 




Figure 1. Part of the tihng Tz/(l/2). 

These tilings all arise from a substitution scheme that is quite similar to the 
pinwheel tiling of Conway and Radin [Rl] . In all cases the tilings have the "sibling 
edge-to-edge" property, which is to say that two daughter tiles of a single parent tile 
can only meet full-edge to full-edge. However, in general the tiling is not globally 
edge-to-edge. Tiles that are not siblings generally do meet in ways that are not 
full-edge to full-edge. 

Chaim Goodman-Strauss [GS] has recently announced that substitution tilings 
with a finite number of prototiles, with each prototile appearing in a finite number 
of sizes, and with the sibling edge-to-edge property, may be associated to local 
matching rules. These rules are less restrictive than local atlasses, and in some cases 
may allow two adjacent tiles to slide freely along their common edge. However, the 
rules are restrictive enough to force any tiling to be hierarchical (i.e., generated by 
the substitution scheme), and hence to be nonperiodic. In a countable number of 
cases considered in this paper, the hypotheses of Goodman-Strauss' theorem are 
met. 



This is striking, because these examples are definitely not globally edge-to- 
edge. Indeed, some of these tilings contain infinitely many pairs of tiles with edges 
that partially overlap. This phenomenon, called "slippage" , can be seen in Figure 
1, especially along the long diagonal. The lengths of the overlaps are not the same 
for all pairs of tiles. In the tiling of Figure 1, in fact, overlaps occur with an infinite 
number of distinct lengths. Thus we have the (seemingly) paradoxical situation 
that Goodman-Strauss' local matching rules force a hierarchical pattern that, in 
turn, forces an infinite variety of local behavior. 

The tilings of this paper also suggest a relaxation of the rules of the tiling 
dynamical systems game. Prior to Conway and Radin's work, a tiling dynamical 
system required that all tiles be generated from a finite set of prototiles by trans- 
lation only, with rotation and refiection not allowed. By those rules, the pinwheel 
tiling, which involves an infinite number of orientations of a single triangle, uni- 
formly distributed about the circle, was not a tiling! The pinwheel example helped 
force a reconsideration of the rules, and a new consensus has emerged that tiles may 
be generated from prototiles by Euclidean motions, not just by translations. 

The present examples suggest a further extension from the Euclidean group to 
the conformal group. In almost all cases, these tilings consist of a single prototile 
appearing in an infinite number of sizes and an infinite number of orientations. 
Although the distribution of sizes is not constant (which it could not be, as the 
conformal group has no Haar measure), it is described by a piecewise-constant 
function. 

Such symmetry is related to properties of the spectrum. In the case of the 
pinwheel, Radin [R2] showed that the statistical rotational invariance of the system 
results in a spectrum that is rotationally invariant, and hence contains no discrete 
component. However, Radin's argument did not address the dependence of the 
spectrum on the radial variable, which corresponds to a length scale. Almost all 
the examples in this paper have a joint distribution of sizes and orientation that 
is rotationally invariant, proving (by Radin's argument) that they have continuous 
spectrum. In addition, in almost all these examples the distribution is absolutely 
continuous in the size parameter. We conjecture that, in these cases, the spectrum 
is purely absolutely continuous. 

Purists who object to the extension to the conformal group may restrict their 
attention to the countable and dense set of examples for which only a finite number 
of sizes appear. They may then take these sizes of the basic triangle as their pro- 
totiles, which then generate all tiles via Euclidean motions. We call such examples 
"rational tilings" . 

The division of the paper is as follows. In §2 we explain the substitution rules 
and the resulting constructions of the tilings, and classify right triangles according 
to whether they generate finite numbers of sizes and orientations, finite sizes and 
infinite orientations, infinite sizes and finite orientations, or infinite sizes and infinite 
orientations. In §3 we examine in detail an example in which tiles appear in two 
sizes and an infinite number of orientations. We prove that this tiling is sibling 



edge-to-edge but exhibits slippage; this tihng contains pairs of tiles that meet in an 
infinite number of distinct ways. In §4 we compute the limiting distribution of sizes 
and orientations in all rational tilings. In §5 we compute the statistical distribution 
of sizes and orientations in irrational tilings. In §6 we examine two exceptional 
rational tilings. 

§2. Definitions and constructions 

As with any substitution scheme, we obtain a tiling of the plane by a succession 
of subdivisions ("defiations"), rescalings and repositionings. We start with a single 
prototile T, define a subdivision rule, and let T^ be the result of subdividing the 
basic tile n times. We will construct infinite sequences of integers A'^i < N2 < ■ ■ ■, 
rescalings Si and Euclidean motions e^ such that ei{si(TN^)) is a proper subtiling of 
ej-(-i(si+i(T/v^^J). By taking the union of the finite tilings ei(si(T/vJ), we obtain a 
tiling of an infinite region, typically the entire plane. 

Our task is complicated by the fact that there is not a straightforward connec- 
tion between A^^ and Si. (Unlike, say, the pinwheel tiling, where Si is an expansion 
by a linear factor 5^"'^.) We cannot, in general, define a simple rule for defiating 
and rescaling by a fixed factor, and simply iterate this rule. We must be more 
subtle. However, in the end, everything does work out. 

Our basic prototile is a right triangle, with base 6, altitude a and hypotenuse c, 
as shown in Figure 2. We refer to the base as the "long" leg, although in principle 
we might have b < a. At this point we put no restrictions on the acute angle 
9 = sin~ (a/c). We divide T into five triangles as in Figure 3. Subtriangles 1-4 
are similar to T, but are smaller by a linear factor of b/2c. Subtriangle 5 is also 
similar to T, but is smaller than T by a linear factor of a/c. This completes the 
first subdivision of T. In the n-th subdivision of T, we take the largest triangles in 
Tn-i, and subdivide each of them into five similar triangles. 




Figure 2. The basic triangle. 

The pattern of subdivision depends on the angle 9. For example, if a > 6/2, 
then the largest triangle in Ti has hypotenuse a, and only this triangle is subdi- 



vided in step 2. Thus in T2 there are 9 triangles, one with hypotenuse a^ jc^ 4 with 
hypotenuse ah/2c, and 4 with hypotenuse 6/2, as shown in Figure 4. The third sub- 
division depends on whether 6/2 > 0? /c, in which case the 4 triangles of hypotenuse 
6/2 are subdivided next, or whether a? /c > 6/2, in which case the one triangle of 
hypotenuse a'^/c is subdivided once again. 




Figure 3. The substitution rule. 




Figure 4. Two subdivisions of a triangle with a > 6/2. 




Figure 5. Two subdivisions of a triangle with a < 6/2. 

The subdivision of a triangle with 6/2 > a is shown in Figure 5. T2 has 21 
triangles, 16 of hypotenuse 6^/4c, 4 of hypotenuse ah /2c, and one of hypotenuse a. 



Finally, if 6 = 2a, then all five triangles in Ti have the same size. All five are 
subdivided in the next stage, yielding T2 with 25 congruent triangles, all of which 
are then subdivided to give T3 with 125 congruent triangles, and so on. This is the 
pinwheel tiling of Conway and Radin [Rl] . 

We will show that, for any angle 9, this subdivision scheme generates nonperi- 
odic tilings of the plane. We first need three technical lemmas: 

Lemma 1: In Tn, the ratio of the hypotenuse of the largest tile to the hypotenuse 
of the smallest tile is strictly less than max(c/a, 2c/b) . 

Proof: The proof is by induction on n. It is clearly true for Ti. Now suppose that 
in Tfc the largest tile has hypotenuse H, while the smallest tile has hypotenuse /i, 
with H/h < max(c/a, 2c/6). After subdividing all the triangles with hypotenuse 
H, the largest hypotenuse is strictly less than H, while the smallest hypotenuse is 
the smallest of /i, aH/c, and bH/2c. The ratio of largest to smallest hypotenuses in 
Tfc-i-i is therefore strictly less than the largest of H/h, c/a and 2c/b. | 

We will often consider an individual tile in Tn , and examine the effect of further 
subdivision of T^ on that tile. When we write t G T„, we mean that t is a single 
tile. When we write t C T^, for some N > n, we mean the collection of tiles in T^ 
whose union is t. When we write that t C T/v is similar to T/j, we mean that the 
individual tiles in this collection fit together to form the single tile t G T^ in exactly 
the same pattern that the individual tiles of Tk fit together to form T. 

Lemma 2: Given an integer n and a tile t E T^, there exists an integer N > n 
such that t C T/v is the union of more than one tile. That is, every tile in Tn, no 
matter how small, is eventually subdivided further. 

Proof: The number of tiles grows without bound (increasing by at least 4 every 
turn). By Lemma 1, the area of each tile is at least min(a^/c^, 6^/4c^) times the 
area of the largest tile. Thus, as A^ ^ 00, the area of the largest tile in T/v must 
go to zero. For A^ large enough, this maximal area will be less than the area of t, 
indicating that the tile t has necessarily been subdivided. | 

Lemma 3: Given an integer n, a tile t E Tn and an integer m, there exists an 
integer N > n such that t C T/v is similar to Tm . 

Proof: The proof is by induction on m. First we consider m = 1. Let A^ be the 
smallest integer such that t G T/v is subdivided. By Lemma 2, such an A^ exists. 
Since t C T/v_i is a single tile, t C T/v is the result of taking a single triangle and 
applying the defiation rule of Figure 3 once. Thus t G T/v is similar to Ti . 

Now suppose the theorem has been proved for m = k, and suppose that t C T/v^ 
is similar to Tk. We consider what happens at the A'o + 1st step. The largest 
triangles in T get subdivided. This means that either the largest triangles in t get 
subdivided, or that none of the triangles in t get subdivided. In the first case t will 
then become similar to T^+i, and we are done. In the second case we consider the 
possibilities for the A'o + 2nd step, and so on. By Lemma 2, the largest triangle in 
t must eventually be subdivided, say on the A^-th turn, so t C T/v will be similar to 
Tjt+i. I 
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Definition: A supertile of order n is a collection of tiles that is similar to T„. 
Equivalently, a supertile of order n is a region of the form e(s(T^)), where s is a 
rescaling and e is a Euclidean motion. 

Theorem 1: Given any right triangle T , there exist tilings of the plane with right 
triangles similar to T , such that any finite set of tiles lies in a supertile, and such 
that the areas of tiles are hounded both above and below. 

Proof: First pick a succession of integers ni, n2, . . . and tiles ti G T^^. Let A^i = ni. 
Pick additional integers A^^, z = 2, 3, . . ., such that, taking tj G Tn^ and subdividing 
T„^ an additional A^^ — rii times, ti C Tj^j. is similar to T^.^-^. By Lemma 3, such 
integers always exist. 

Now pick a triangle similar to Tq and place it in the plane. One can construct 
a supertile Si of order A^i that contains this triangle in position ti. One then 
constructs a supertile 5*2 of order A^2 such that Si sits inside 5*2 as t2 sits inside 
T„2. One continues the process, building supertile Sk+i such that 5'^ sits inside 
Sk+i as tk+i sits inside T^^+i. 

The union of all the supertiles is a tiling of an infinite region. For almost all 
choices of the t^'s (e.g., having the edges of ti lie in the interior of T^j,. infinitely 
often), this region will be the entire plane. 

Since the ratio of largest to smallest triangle is uniformly bounded for 5'^^ 
by Lemma 1, no tile may have hypotenuse longer than max(c/a, 2c/6) times the 
hypotenuse of a fixed tile, and no tile may have hypotenuse less than min(a/c, b/2c) 
times the hypotenuse of the same fixed tile. This provides both an upper and lower 
bound to the size of the tiles. | 

In this construction, many choices were made. Different choices generally lead 
to different tilings, but these different tilings have many properties in common. 
The following Theorems 2 and 3 apply to all tilings constructed in the manner of 
Theorem 1. 

Theorem 2: In a tiling, the number of orientations in which the basic triangle 
appears is finite if Q Jtx is rational and infinite if 6 /ti is irrational. 

Proof: First suppose that 9 /ti is rational. Consider a tile, positioned as in Figure 
3, with side h along the x axis. Let P denote refiection about the x axis, and let 
Ra denote a counterclockwise rotation by angle a. After subdividing once, the 
orientations of the five daughter tiles, relative to the parent tile, are given by the 
following elements of 0(2): RgP (twice), Rg, Rn+e, and RT^/2+eP- The orientations 
of tiles in a further subdivision are words in these five elements of 0(2). However, 
with O/tt rational, these five elements generate a finite subgroup of 0(2), so only a 
finite number of orientations can ever appear in a future subdivision. 

Since any region of our tiling of the plane sits inside a supertile, any two tiles 
must have their orientations, relative to the supertile itself, in this group. Thus 
their orientations, relative to each other, and hence to a fixed reference tile, must 
lie in the group. Thus only a finite number of orientations can appear in the tiling. 



Now suppose that 6 /ti is irrational. In a basic subdivision, we will keep track 
only of the 4 triangles of hypotenuse 6/2, ignoring the triangle of hypotenuse a. 
When these 4 triangles divide, we will only keep track of the 16 resultant triangles 
of side 6^/4c, and so on. In the second generation we find orientations 1 = [RgP)"^ 
and -R29 = -R^' among others. In the 2n-th generation we find 1, R2e, Rab, . . . , R2ne- 
Since 9/n is irrational, these 2n + 1 orientations are distinct. Since our tiling of the 
plane contains supertiles of arbitrarily large size, there is no bound to the number 
of different orientations that appear. | 

Theorem 3: In a tiling, the number of sizes in which the basic triangle appears 
is infinite i/ln(sin(6'))/ ln[cos(6')/2] is irrational and finite z/ln(sin(6'))/ln[cos(6')/2] 
is rational. In particular, z/ ln(sin(6'))/ ln[cos(6')/2] = p/q, with p and q relatively 
prime integers, then the number of sizes in the tiling is max(p, q) . 

Proof: Let A = a/c and let B = b/2c. If ln{A)/ln{B) is irrational, the only way 
two monomials A°^B^ and A'^B'^ can equal is if a = c and b = d. We will show 
that the sizes of triangles in T^ (relative to the original triangle) is given by such 
monomials, and that the number of distinct powers of A grows without bounds as 
n ^ 00. This will show that the number of distinct sizes grows without bound a 
n ^ 00. 

In each subdivision there are 4 tiles of size B relative to the parent and one 
tile of size A. Thus the descendants of a given tile all have sizes that are monomials 
j[aj^b j-gia^ive to the ancestor. For every n > 0, T^ contains at least one tile with 
size A^B^; just take a B child of a S child of ... of one of the original B children 
(or the B child itself, if it has not subdivided). For every n > 0, T^ contains at 
least one tile with size A^B^; take a BBB . . . descendant of the original A child. 
Once n is large enough to have the original A child divide, there is at least one 
tile with size A^B^. In general, once n is large enough to allow a a-th generation 
AA . . . child, it will always have at least one tile whose size has exactly a powers of 
A. This completes the irrational case. 

Now suppose that ln{A)/ln{B) = p/q, with p and q relatively prime. Thus 
j^q — j^v Lg^ J- = A^'P = B^''^. Every monomial A'^B^ is a power of r. Assume for 
the moment that p < q. By Lemma 1, the ratio of sizes of any two tiles is greater 
than B = r'^. Thus only at most q distinct sizes can appear. To see that q sizes do 
appear, we note that A^S^'o, A^B^^ , . . . , A'^~^ B^''-^ are all distinct powers of r. 

li p > q, Lemma 1 states that the ratio of any two sizes is at least A = r^, 
so at most p different sizes can occur. We produce p different sizes by examining 
different powers of B. In either case, the number of distinct sizes is max(p, q). | 

We refer to tilings with ln{A)/ ln{B) = p/q as (p/q) rational tilings, and denote 
the class of such tilings as Til{p/q). For z irrational, we will similarly denote the 
class of tilings with ln{A)/ ln{B) = z as Til(z). The different tilings in a class are all 
derived from the same substitution rule, and have many properties derivable from 
this rule. When discussing such properties, we will sometimes refer to a typical 
element of the class as "the tiling Til{z)" . 



Note that ln(sin(6'))/ln[cos(6')/2] is a strictly decreasing function of 6 on the 
interval (0, 7r/2). From this monotonicity, and from the countability of the rationals, 
it is clear that only a countable set of angles 9 give rise to a finite number of 
rotations, and only a countable set of angles 6 give rise to a finite number of sizes. 
The intersection of these two countable sets turns out to be a single point. 

Theorem 4: The only angle that gives rise to both a finite number of orientations 
and a finite number of sizes is 9 = it/ 4. That is, the tiling Tz/(l/3). 

Proof: Let x = exp(i6'), with < 9 < 7r/2. We are looking for solutions to the 
equation sin(^)'^ = [cos(^)/2]^, which we rewrite as 

2'i{x + xf = 2^P{-iy{x-xy, (2.1) 

where x = exp{~i9) = x~^. Note that, for fixed p, q, there is at most one solution 
to equation (1) in the first quadrant, since ln[sin(6')]/ ln[cos(6')/2] is monotonic. 

Since by assumption ^ is a rational multiple of tt, x is a primitive n-th root of 
unity for some integer n. If x is a solution and q is even, then equation (1) has (real) 
integer coefficients, and all the primitive n-th roots of unity are also solutions. If q 
is odd, all the primitive n-th roots of unity are solutions either to equation (1) or 
to the conjugate equation 

2'^{x + xf = 2^Pi'^{x-xy. (2.2) 

Equation (2), with q odd, has no solutions in the first quadrant, as the right hand 
side is positive but the left hand side is negative. Since equation (1) admits only 
one solution in the first quadrant, there must be exactly one primitive n-th root of 
unity in the first quadrant. 

This means that n must equal 5, 6, 7, 8, 10, 12 or 18. Checking these 
individually, we see that only n = 8, or ^ = 7r/4, yields a rational value of 
ln[sin(6')]/ln[cos(6')/2]. In that one case sin(6') = cos(6') = ^2/2, and sin^(6') = 
cos(^)/2. I 

In a periodic tiling, all of the sizes and orientations are exhibited in a compact 
region, so neither the number of sizes nor the number of orientations can be infinite. 
Thus we have 

Corollary If z ^ 1/3, then the tiling Til{z) is not periodic. 

In fact, it will turn out that Ti/(l/3) is not periodic, either. This will be shown 
in §6. 

§3. An example with two sizes. 

In this section we consider in detail the tiling Ti/(l/2), shown in Figure 1. 
This example is chosen not as a special case, but rather as a simple example of 
some general phenomena. Based on the statistical analysis of §4, we expect all 
tilings Til{p/q), with q > 1 and p/q ^ 1/3, to be qualitatively similar to Til{l/2). 
Specifically, in all these cases the population matrix has two or more eigenvalues 

9 



with modulus bigger than one. This causes fluctuations in the statistical composi- 
tion of supertiles Sn to grow with n. This, in turn, can cause the tiling to fail to be 
globally edge-to-edge. 
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Theorem 5: Til{l/2) is a tiling with two sizes of tiles, each of which appears in 
an infinite number of orientations. The substitution scheme has the sibling edge-to- 
edge property, but the tiling is not globally edge-to-edge. Specifically, the tiles meet 
in an infinite number of ways. 

Proof: The existence of the tiling, the number of sizes and the number of orientations 
follow from Theorems 1-3. The sibling edge-to-edge property is manifest, if we 
consider the basic tile to have four vertices - the three obvious ones and the midpoint 
of the long leg. The difficulty is in proving that tiles meet in an infinite number of 
distinct ways. 




Figure 6b 




Figure 6c 



Figure 6. Three stages of subdivision for Til{l/2). 

The process by which this happens is as follows. There are certain special lines 
in the tiling. The long diagonal in Figure 1 is an example. On each side of such a 
line there are triangles, all of the same size, whose hypotenuses or long legs make 
up part of the long line. The pattern is different on the two sides of the line, with 
one side having (say) more legs in a certain region and the other side having more 
hypotenuses. This imbalance causes the tiles on one side of the line to appear shifted 
relative to those on the other side, a phenomenon we call "slippage" . We will exhibit 
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regions where the imbalance is arbitrarily large. Since the imbalance is unbounded, 
the slippage of one side relative to the other reaches arbitrarily high multiples of 6, 
modulo c, where b and c are the lengths of the long leg and hypotenuse. Since b/c 
is irrational, this gives rise to an infinite number of ways in which one triangle can 
meet another across the long line. 

Lines where slippage occurs are called "fault lines" . Note that each fault line 
has only finite length, and allows only a finite amount of slippage. However, we 
will find fault lines of arbitrarily long length with arbitrarily much slippage. This 
precludes there being a finite bound on how many ways one triangle can meet 
another. 

The Tz/(l/2) tiling is based on a right triangle with legs a — W2(vT7 — 1) and 

b = y/Vf — 1 and hypotenuse c = 4. Several degrees of subdivision are shown in 
Figure 6. At each level there are two sizes of triangle, whose linear sizes differ by 
a factor of a/c. The heavily shaded lines in Figure 6 are fault lines. For each T^, 
with n even, only large triangles abut the illustrated fault line. Once a fault line is 
formed, the triangles on opposite sides of the line evolve separately, and begin to 
slip. We shall prove that this slippage increases without bound. 

Ti/(l/2) may be viewed as a traditional substitution system, with two pro- 
totiles, which we call B and S (for big and small). Each subdivision, followed by 
rescaling by c/a, may be viewed as a replacement of each S triangle by a S triangle, 
and replacement of each B triangle by a S triangle and four 5" triangles. Let the 

population of T„ be \1/^ ~ ( /\r ) ' where A''^ and Ns are the numbers of big and 

small tiles. \1/ satisfies \I/n+i = Af ^n? where the population matrix is 

M={1 I), (3.1) 

with eigenvalues X± = (l±\/T7)/2 and eigenvectors (± = I 1 . As n grows. 



the ratio of A^ to Ns approaches (vl7 + l)/8 ~ 0.6404. The exact populations are 

(3.2) 

Ns{n)=4NLin-l) 



Note that |A_| > 1, so that |8Ab - (Vl7+ l)As| grows with n. 

Next we consider what happens along a fault line. To do this we must consider 
the boundary of T2n- Note that the hypotenuse and long leg of T2 consist only of 
hypotenuses and long legs of big triangles. Applying the subdivision again, we get 
that the hypotenuse and long leg of T4 also consists only of hypotenuses and long 
legs of big triangles. Similarly for all T2n- 
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The evolution of these legs and hypotenuses is a substitution system in its 
own right, only in one dimension. There are four symbols, H^ , H~ , L+ and L~ , 
representing the two orientations of the hypotenuse and long leg, respectively. From 
Figure 6b, we see that the substitution rule, which we denote ao, is 

ao{H+) = L+ L- H+ , aoiH') = H' L+L' , ao{L^)=H^H^. (3.3) 

Since L"*" and L~ only appear in the combination L~^L~ , we can define a new symbol 
L = L^L~ and have a substitution system with 3 elements, whose rule we denote 
a: 

a{H+)=LH+, a{H-)=H-L, a{L) = H+H+H-H' . (3.4) 

Lemma 4: The sequence a'^{H^) contains neither the subsequence LL nor the 
subsequence H~H~^ . 

Proof: The proof is simple induction. The only way to generate an LL is from an 
H~H^ , and the only way to generate an H~ H^ is from LL. Since neither appear 
in the first generation, neither appears in any subsequent generation. | 

Lemma 5: The sequence a'^{H~^) does not contain a subsequence of more than 6 
consecutive H 's. 

Proof: Since LL does not occur in a^~^{H~^), the longest possible sequence of if's 
in a'^{H^) would come from a sequence H~^LH~ in a'^~^{H~^). This gives rise to 
LH+H+H+H-H-R-L, or 6 H's in a row. | 

Let f{n) equal the number of L's in the first half of the sequence a'^{H~^) minus 
the number of L's in the second half of the sequence. As we shall see, f{n) is closely 
related to the extent to which slippage occurs along the largest fault line in T2n+2- 

Lemma 6: // |/(n)| > 6, then \f{n + 1)| > |/(n)| + 2. 

Let s and s' denote the first and second halves ofa'^{H^), respectively. Suppose 
that s contains h H^s and / L's, while s' contains h' H^s and /' L's. Note that 
/(n) = I — I' = h' ~ h. Since each H generates an H and an L, while each L 
generates four LT's, a{s) contains 41 + h if' s and h L's, while cr(s') contains 41' + h' 
if s and h' L's. Thus a{s) contains 2f{n) more terms, but /(n) fewer L's, than 
a{s'). 

Now suppose f{n) > 0. The first half of a'^~^^{H~^) is all of cris), minus the 
last f{n) elements, while the second half of a'^^^{H^) is the last /(n) elements of 
a{s) and all of cr(s'). Thus /(n + 1) equals —fin) minus twice the number of L's 
in the last f{n) elements of (t{s). Since f{n) > 6, there must be at least one L in 
the last /(n) elements of o"(s), so \f{n + 1)| = —f{n + 1) > f{n) + 2. 

If f{n) < 0, then the first half of a'^+^{H+) is all of (t(s), plus the first |/(n)| 
elements of cis'). We then have /(n + 1) equalling —f{n) plus twice the number 
of L's in the first |/(n)| elements of (j{s). By Lemma 5, there must be at least one 
such L, so \f{n + 1) = /(n + 1) > 2 - /(n) = |/(n)| + 2. ■ 

Lemma 7: lim^^oo |/(^)| = +oo. 
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Proof: By explicit computation, /(I) = 1, /(2) = -1, /(3) = 1, /(4) = -3, 
/(5) = 3, /(6) = -5 and /(7) = 9. By Lemma 6, for n > 7, |/(?i+ 1)| > |/(n)|, so 
|/(n)| > n + 2 goes to infinity as n ^ oo. | 

We have proven that |/(?i)| grows without bound, which is aU that we need. In 
fact, |/(n)| grows exponentially. For large n, an approximate fraction 2/(3 + VT?) 
of the elements of a''{H+) are L's, so /(n + 1) ^ -f{n) + 2/(n) x 2/(3 + v^) = 
/(n)(l — \/T7)/2. The growth rate, (1 — vT7)/2, equals A_, the second eigenvalue 
of the population matrix M. 

We now return to the question of slippage along fault lines. Consider two large 
triangles that meet hypotenuse to hypotenuse to form a rectangle, as in Figure 7. 
Let P and R be the ends of the common hypotenuse, and let Q be the midpoint. 
Rotation by tt about Q sends each triangle into the other. 



Triangle 1 




Figure 7. A fault line. 

Subdivide the pair of triangles 2n times. The pattern of subdivision of triangle 
2 along the interval PQ is the same as that of triangle 1 along the interval RQ. We 
have already seen that only large triangles abut the main diagonal, and that they 
do so only along their hypotenuses and long legs. For each point x on the interval 
PQ, let gn{x) be the number of complete long legs on the triangle 1 side of Px 
minus the number of complete long legs on the triangle 2 side of Px. 

Lemma 8: \gn{Q)\ > f 1/(^)1 ~ 1- -^'^ particular, Limn^oo\9n{Q)\ = oo. 
Proof: Consider the the 2n-fold subdivision of the pair of triangles in Figure 7. The 
hypotenuse of triangle 2 gets divided into a^{H~^). Let M„ be the point on the 
hypotenuse corresponding to the middle of the sequence a'^{H^). By construction, 
there are 2f{n) more long legs, and f{n) fewer hypotenuses, between P and M^ 
than between M„ and R. Since legs have length 6, hypotenuses have length c, 
and 26 > c, M„ lies a distance f{n){2b — c) closer to R than to P, or a distance 
/(n)(26-c)/2 beyond Q. 

Suppose /(n) < 0. Then Mn lies between P and Q. The number of hypotenuses 
between P and Q is at least the number of hypotenuses between P and M„, and so 
is at least —f{n) more than the number of hypotenuses between Q and R. Since the 
length of PQ equals that of QR, and length equals bx legs plus ex hypotenuses, 
there are at least —{c/b)f{n) more legs in QR than in PQ. Thus \gn{Q)\ is at least 
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the integer part of (c/6)|/(n)|, which is greater than {c/b)f{n) — 1. 

If f{n) > this situation is reversed. Then M^ hes between Q and R, and the 
interval between M^ and R has the surplus of f{n) hypotenuses. The interval QR 
has a surplus of at least that many hypotenuses, so the interval PQ has a surplus 
of at least {c/b)f{n) legs. Taking the integer part, we see that PQ has a surplus of 
at least {c/b)f{n) — 1 complete legs. | 

We are now in a postition to prove Theorem 5. Let a; be a point on PQ that is 
the vertex of a tile in the subdivision of triangle 2, and let y be the nearest vertex, 
on PR, of a tile in the subdivision of triangle 1. The interval xy is the interval of 
contact between two tiles, one in triangle 1 and the other in triangle 2. If y lies 
between P and x, then the length of xy is congruent, modulo c, to bgn{x). If x lies 
between P and y, then the length is congruent either to —bgn{x) or b — bgn{x) (mod 
c) , depending on whether xy is part of the hypotenuse or long leg of a tile on side 
1. 

As we pick points x from P to Q, gn{x) goes, by steps of one, from to gn{Q)- 
Thus the intervals xy take on at least \gn{Q)\/2 distinct lengths. Since \gn{Q)\ 
grows with n, all we need for Theorem 5 is to show that fault lines modeled on 
a'^{H~^), for arbitrarily large values of n, occur in the tiling Ti/(l/2). 

By construction, Til{l/2) contains supertiles modeled on T^, for arbitrarily 
large values of m. If m is even, then the primary fault line exhibited in Figure 6 is 
modeled on a^{H~^), with m = 2n + 2. If m is odd, then there is also a supertile 
modeled on T^-i, namely the descendants of the large tile in the first subdivision 
of T. Since ?7i — 1 is even, this supertile contains a fault line modeled on a^{H~^), 
with m = 2n + 3. Since m is unbounded, we have obtained our requisite arbitrarily 
long fault lines. 

This completes the proof of Theorem 5. | 

§4. Statistical properties of rational tilings 

In this section we consider the statistical distribution of sizes and orientations 
of tiles in the various rational tilings Til{p/q). More precisely, we consider the 
distribution of sizes and orientations in a supertile modeled on Tn, and take the 
limit as 71 ^ oo. We first show that the distribution of sizes approaches a simple 
limit as n ^ oo. The limiting distribution is given by the eigenvector associated 
to the largest eigenvalue of a, p x p (or q x q) population matrix. We also analyze 
the second eigenvalue of this matrix. The failure of the edge-to-edge property for 
Til{l/2) was a result of fiuctuations that were governed by this second eigenvalue. 
We conjecture that the edge-to-edge property holds only for those tilings with second 
eigenvalue smaller than 1, and we classify these tilings. 

We then turn to the joint distribution of sizes and orientations in rational p/q 
tilings. We show that, for each size, the distribution of orientations is asymptotically 
uniform. Specifically, we parametrize 0(2), the group of orientations, by two copies 
of the unit circle. Given an interval in this set, the fraction of tiles in T^, of a 
given size, whose orientations lie in that interval, approaches a constant times the 
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length of the interval. In the terminology of Radin [R3], the tilings T{p/q) exhibit 
"statistical rotational symmetry" . 

By Theorem 3, Til{p/q) contains triangles of m = max(p, q) distinct sizes. 

Til{p/q) is equivalent to a traditional substitution tiling, with prototiles Di, . . . Dm 

of TO sizes. We take Di to be the largest size and D^ to be the smallest. Subdivision 

and linear rescaling by r~^, where r = {a/cY''P = (6/2c)^''^, takes -D^+i to Di, and 

takes Di to four copies of Dq and one copy of Dp. That is, the population matrix, 

which gives the population of T^+i in terms of the population of T„, has matrix 

elements 

M ifj = z + l; 

^1 1 ifj = 1 and i=p; ^^_^^ 

14 if J = 1 and i = q; 

V otherwise. 

The properties of M are summarized in the following theorem: 
Theorem 6 

1) The characteristic polynomial of M is 

v(X) -l^'- ^'"' - 4 '^P< ^' (4 2) 

P^-^)-\xP-4XP-'i-l tfp>q. ^ ' 

2) The largest eigenvalue of M is r~^. 

3) There are exactly q eigenvalues with modulus greater than one. 

4) The eigenvectors i/j of M , for fixed eigenvalue \, take the form 

^^ = \k - \k-PH{k -p-l)- 4X''-iH{k - Q - 1), (4.3) 

where H{n) is the discrete Heavy side function 

H{n) = {^ '^''-°' (4.4) 

1 otherwise. 

5) Asymptotically, the number of tiles of size Dk is a fraction 

Vk - ^-^{a'Hip -k) + h'Hiq - k))r-^'' (4.5) 

4c^ 

of the total. 

6) Asymptotically, the area covered by tiles of size D^ is a fraction 

a^H{p-k)+9H{q-k) 

Pk = 2 I h2 ^-6 

of the total. 

Note that, if p < g, then H{k — g — 1) is identically zero and H{q — k) is 
identically one. If p > q, then H{k — p — 1) is identically zero and H{p — k) is 
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identically one. As written, expressions (4.3), (4.5) and (4.6) apply to both the 
p < q and p > q cases. 

Corollary: Let t be a tile of size k. Let Nk{n) be the number of descendants oft 
after the n-th application of the substitution rule. Then 

Proof of Corollary: From the distribution (4.5), we compute the average area per 
tile to be (1 — r'^)ah{j>a^ + qb'^)/8c'^. A tile of size D^, subdivided and rescaled n 
times, has area a6r^^~^"/2. Dividing by the area per unit tile we obtain (4.7). | 

Remark: If q > 1, then the second largest eigenvalue of M is greater than one. 
The fluctuations in population associated to the corresponding eigenvector then 
grow with subdivision, although they do not grow as fast as the population as a 
whole. It was precisely this phenomenon that caused the edge-to-edge property to 
fail for Til{l/2). If g = 1, then the second eigenvalue is less than one. In §6 we 
shall see that, in Til (2), this causes the triangles to meet in only a finite number of 
ways. 

Conjecture: If p and q are relatively prime integers, q > 1 and p/q ^ 1/3, then 
the tiling Til{j>/q) has tiles that meet in an infinite number of different ways. If 
p> 1, then the tiles in Til{j>) meet in only a finite number of ways. 

Proof of Theorem 6: For k < m, the k-th row of the vector equation Mt/j = Xi(; 
reads 

•0fc+i = Xi^k - {Sk,p + 4:dk,q)i^i. (4.8) 

Setting ifji = A, we repeatedly use (4.8) to obtain expression (4.3) for 'i/'2, • • • V'm- 
Plugging this into the m-th row of Mt(; = Xtp then gives the characteristic polyno- 
mial (4.2). 

Now suppose p < q, and consider the function p{\) for A real and positive. 
Note that p{r-'^) = r''^'^ - r'^P-'^i - 4 = Ac'^/b'^ - Aa^ /b'^ -4 = 0, since a'^ + b'^ = c^. 
When A > 1, p'{\) = qX'^'^ - {q - p)Xi-P-^ > 0. Thus p{X) > for all A > r'^. 

Now let A 7^ r~^ be a root oi p{X). We will show that |A| < r~^. If A is real 
and positive, then A < r~^. If A is not real and positive, then, since q and q —p are 
relatively prime, we cannot have A*^ and A^~^ both real and positive. Thus, by the 
triangle inequality, 

= \p{X)\ = |A« - X1-P - 4| > |A'?| - |A«-P| - 4 = p{\X\). (4.9) 

Since p(|A|) < 0, we must have |A| < r~^. 

We count the number of roots in the unit circle via the argument principle, 
tracking the argument of p{ex.p{i9)) as 9 goes from to 2%. When |A| = 1, p{X) 
always has negative real part, as 4 > jA"^ + X'^~p\. Thus the winding is zero, and 
none of the roots of p{X) lie in the unit circle. Thus there are q roots, counted with 
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multiplicity, outside the unit circle. Since p'(A) is never zero outside the unit circle, 
all these roots are distinct. 

Now suppose p > q. Then p{X) = A^~^(A^ — A^~^ — 4). By the same arguments 
as before, p(r~^) = 0, and p{X) > for real A > r~^. For any eigenvalue A other 
than r~^, = \p{X)\ > p{\X\). Since p(|A|) < 0, |A| < r~^. By the argument 
principle, there are p — q roots inside the unit circle, since the dominant term of 
p{X) on the unit circle is 4X^~'^. This leaves q roots outside the unit circle. This 
completes the proof of statements 1-4. 

To obtain the asymptotic distribution of sizes, we must decompose the ini- 
tial population into eigenvectors of M. The asymptotic distribution will be the 
eigenvector corresponding to the largest eigenvalue, assuming the coefficient of that 
eigenvector is nonzero. We have shown that this largest eigenvalue is r~^. Since the 
total area of the system grows by a factor of r~^ each time, the coefficient of this 
eigenvector is not zero. Thus the asymptotic distribution of population is given by 
a multiple of the expression (4.3), with A = r~^. Normalizing, we obtain expression 
(4.6). Multiplying this by the area of each tile and normalizing again gives the 
asymptotic distribution of areas (4.7). | 

We now turn to the joint distribution of sizes and orientation. To do this we 
must first parametrize the space of possible orientations of a single size. This space 
is isomorphic to 2 copies of the unit circle. We specify both the handedness of the 
triangle and the direction a fixed vector in it points in the plane. We take as our 
reference vector the ray from the small angle to the right angle. In Figure 8, the first 
triangle has orientation (+,(j)i), while the second has orientation (— ,(/)2). We will 
let the O denote the ordered pair (±, (/>), and let dQ = d(/)/47r be the Haar measure 
on the space of orientations. The space of all possible tiles up to translation, which 
we denote X, is 2m copies of 5"^. 



Figure 8. Orientation of triangles. 

Let / be a function on X . Given a collection S of tiles, let < /, 5" > be the 
average value of / on the individual tiles in S, where each tile is given equal weight. 
Let < f, S >' be the average value of / on the individual tiles in S, where each tile 
is weighted proportionally to its area. 

Theorem 7 Assume a tiling Til{p/q), with p and q relatively prime andp/q ^ 1/3. 
Let f be a continuous function on X , and let {Sn} he a sequence of supertiles of 
increasing size. Let dv = J2k ^kdVtk and dp = J2k Pkd^k be measures on X , where 
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z/fc and pk are as in Theorem 6. Then 

lim < /, Sr, >= I fdu (4.10a) 



X 

lim <f,Sr,>'= [ fdp. (4.106) 

These limits also apply if f is the characteristic function of an interval in X . 

Proof: We first reduce the problem to establishing (4.10a) for an arbitrary con- 
tinuous function /. Once we have established (4.10a), (4.10b) follows by applying 
(4.10a) to the continuous function /(/c,0) = r^^/(A;,0). Once equations (4.10ab) 
have been established for continuous functions, the extension to characteristic func- 
tions is standard (for details see [CFS]). Let / be an interval, and let xi be its char- 
acteristic function. We choose continuous functions f^, such that /~ < X/ < /J^, 
and such that Xivn^^Q j.^ ff^du = Jjdv. lim^^cx) < Xi^Sn > is sandwiched be- 
tween lim^^oo < /e~5 'S'n > and lim^^cx) < /J^, Sn >, hence between f^ fT^^ ^^^ 
Ix ftdi', and so must equal /^ df. 

To establish (4.10a) we must introduce some notation. If / is a function on X 
and 5" is a collection of tiles, let (/, S) be the sum of / evaluated on the individual 
tiles of 5". Let $ denote the action of subdividing and rescaling. That is, $ acting 
on a tile of size -Dfc+i gives a tile of size D^ of the same orientation, while $ acting 
on a tile of size Di gives one tile of size Dp and four of size Dq, having various 
orientations. Let $* be the dual of $ by (■, ■), acting on the space of functions: 

($7,^) = (/,$(^))- (4.11) 

Let fn = ($*)'^/. Note that $* is linear and sends non-negative functions to non- 
negative functions, so if / < g, then /„ < Qn- 

Now suppose we have a sequence of tiles in and supertiles 5"^ = ^^tn- We have 
that 

# of tiles in 5'^ # of tiles in 5"^ 
Since the number of tiles is given asymptotically by (4.7), (4.10a) is equivalent to 



r 



2n 



fn converging uniformly to J^ fdu times 



We examine the spectrum of the linear operator $* on the function space C{X). 
The key lemma, whose proof we defer, is 

Lemma 8 The spectrum of $* is pure point. Co is an eigenfunction with eigen- 
value r~^. All other eigenvalues have norm strictly less than r~^. Any continuous 
function f can be written as a (possibly infinite) sum of eigenf unctions 0/ ($*)l'^~^l, 
such that a subsequence of partial sums converges uniformly to f . 
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Remark: If p < q, it turns out that there are a number of functions ( for which 
(^$*)l9-p|^ = 0, but $*C 7^ 0. Thus, to achieve a basis for the space of continuous 
functions, we must use eigenfunctions of ($*)I9-pI rather than just eigenfunctions 
of $*. 

Given the lemma, we write 

CXD 

/ = $]qO. (4.14) 

Since a subsequence of the partial sums converges uniformly, and since Co has a 
positive minimum, for each e > we can find an integer N such that 

/-= |^c,Oj-<o</< [X]cfcOj+<o = /+, (4.15) 

where each Q is an eigenfunction with eigenvalue A^. (Strictly speaking, Q is merely 
an eigenfunction of ($*)I9~pI, not necessarily of $*, but this distinction makes no 
difference). Since |r~^Ai| < 1 for alH > 0, for n > \p — q\ we have 

N N 

(CO - e)Co + J2ck{r^^,rC^ < ^^'"fn < (CQ + e)Co + J]cfc(r2A,)"C^• (4.16) 

As n ^ oo, the left hand side converges uniformly to (cq — e)Co, while the right hand 
side converges uniformly to (cq + e)Co. Since e is arbitrary, r^"/n must converge 
uniformly to CqCo- 

All that remains is to compute cq in terms of /. Since Co is invariant under 
rotation and refiection, cq must be of the form ^^ dk J f{k, Q)dQ for some universal 
constants dk- By comparing characteristic functions of different sizes, we see that 
the constants dk must be proportional to z/fc. Finally, for the constant function 
/ = 1, fn{k,Q) is the number of descendants of a tile of size Dk, which we have 
already computed in (4.7). This fixes the proportionality constant. | 

Proof of Lemma 8: $* commutes with rotations, so we may simultaneously diago- 
nalize $* and the rotation operator —i-^-l. The eigenvalues of the rotation operator 
are of course the integers, with each eigenspace being 2rn.-dimensional. Specifically, 
the eigenspace corresponding to an integer n is the span of the 2m functions ob- 
tained by restricting exp(m(/)) to each of the 2m circles in X. Operators on finite 
dimensional spaces always have pure point spectra. Summing over n, we see that 
the spectrum of $* is pure point. 

On each 2m dimensional subspace corresponding to the Fourier mode exp(m(/)), 
the action of $* is described by a 2m by 2m matrix E, which we write as an m x m 
array of 2 x 2 matrices. Let 6 = tan~^ (a/c) be the acute angle in our basic triangle. 
Consider the matrices 

/ exp[-m(^ + 7r/2)]\ 

Vexp[m(6' + 7r/2)] J' ^ ' 
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B 



(4.18) 



exp{in9) + exp{in{9 + tv)) 2 exp{—in9) 

2 exp(in6') exp(— in6') + exp(m(— 6* + tt)) 

A and S describe the orientations of the five daughter tiles in terms of the orienta- 
tion of the parent tile, as expressed in the n-th representation of the rotation group 
SO{2). Specifically, A describes the daughter tile of hypotenuse a, while B describes 
the four daughter tiles of hypotenuse 6/2. In each case the first column describes 
the daughters of a positively oriented tile, while the second column describes the 
daughters of a negatively oriented tile. See Figure 9. 





Figure 9. Orientation of daughter tiles. 
Our matrix E has the following matrix elements: 



E. 



«j 



(4.19) 



if j = z- 1; 

if i = 1 and j — p; 

if i = 1 and j — q; 

otherwise. 

Notice that E is essentially the transpose of our population matrix M, with the 
daughters represented by the matrices A and B rather than the numbers 1 and 4. 
We compute the eigenvalues and eigenvectors of E by directly solving the equation 
Etp = >4. 
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Suppose V'm = V, where t; is a 2 component vector. Each row but the first of 
Ei/j = Xijj implies that if^k-i = AV'fc, and hence that i/jk = y^~^v for all k. The first 
row then says that 

(A"^ - W^-'iB - r'-PA)v = 0. (4.20) 

Taking the determinant of the matrix on the left hand side, we obtain the charac- 
teristic polynomial of E, 

p( A) =A2"^ - A^"^-^ (2 cos(n6') + 2 cos(n7r + n6')) + A^'^-^^ (2 cos(n7r) - 2) 

_ ^2m-2p _ ^2m-p-q^^ COs(n7r/2)). ^^'^^^ 

First we consider n = 0. If p < q, p{X) = (A*? - X^'P - 4)(Ap + 1)Xi-p. We 
recognize the first factor as the characteristic polynomial of M, with one root r~^ 
and all other roots smaller in norm. The roots of the second and third factors have 
norm 1 and 0, respectively. The eigenvector corresponding to eigenvalue r~^ is 



•0/0=^2^(^1, i.e., Co (up to scale). 



We can understand the eigenvalues as follows. The roots of the first factor 

all correspond to t; = I J , and describe fiuctuations in the numbers of tiles of 

various sizes, irrespective of orientation. This is the problem we previously studied 

in Theorem 6. The unit and zero eigenvalues correspond to v = i ^ ) . The 

zero eigenvalue has algebraic multiplicity q —p and geometric multiplicity 1. Since 
the basic subdivision produces equal numbers of positively and negatively oriented 
"B" tiles, after q — p subdivisions there must be the same number of positively 
and negatively oriented tiles of each size Dp+i, . . .Dm- Since there is but one 
"A" daughter of each subdivided tile, the imbalance between positive and negative 
orientation in the larger sizes neither grows nor shrinks, but oscillates with period 
2p. These correspond to the roots of A^ + 1. 

If p > q, then p{X) factorizes as p{X) = (A^ - 1 - 4Ap~'^)(Ap + 1). Again, the 
first factor is the characteristic polynomial of M, and governs the total number of 
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tiles of each size, with all eigenvectors having v = i ^ 1 . The largest eigenvalue is 
r~2, with eigenvector Co- The roots of (A^ + 1) describe oscillations in the numbers 
of positively vs. negatively oriented tiles and have v = 



-1^ 

Next we consider n odd. Then exp(m7r) = —1, and p{X) simplifies to A^"^ — 
^^2m-2q _ ^2m-2p xhis is just the characteristic polynomial of M applied to A^. 
By Theorem 6, the largest roots have A^ = r"^, or A = ±r~^. 

Finally we consider n even and nonzero. Thenp(A) = X'^'^ — X'^"^~'^{4cos{n9)) — 
y2m-2p _ x^^-p-Q (^^^QQ^^jm l2)) . We show that all roots are smaller than r"^ by 
the argument principle. Note that 6 is an irrational multiple of tt, so cos{n6) ^ ±1. 
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On the circle |A| = r ^ we have 

|A2"^-^(4cos(n6'))| + \X^'^-^P\ + \X^'^-P-'i{4cos{nn/2))\ 

_i >2m| f b^\cos{ne)\ , a^ ^252 



(4.22) 



C2 



<|A 



2m I 



Thus the A2"^ term of p{X) dominates on the circle |A| = r~'^, so the winding of the 
argument of p{r~'^ ex.p{ia)), as a goes from to 27r, is 2m, and all 2m eigenvalues 
of E lie inside the circle of radius r~^. 

This completes the analysis of the spectrum of $*. Now we need only consider 
the decomposition into eigenvectors. Since / is continuous, it has an absolutely 
summable Fourier series. So we may write /(/c, O) = ^^CnV'fc 6xp(m(/)), where 
the c„'s are absolutely summable and ipjT' is a vector in C2"^ with no component 
larger than 1. In particular, the partial sums converge uniformly to /. But \l/fc 
is itself a sum of eigenvectors of i?''^"^!, so we may rewrite our sum as a sum of 
eigenfunctions of ($*)I9~pI. | 

§5. Irrational tilings 

In this section we consider tilings Til{z), with z irrational. The analysis is 
formally similar to that of §4, except that we are now dealing with an infinite number 
of possible sizes. In place of the discrete size parameter k we introduce a continuous 
size parameter s. In place of the discrete evolution operator ^^ we introduce a 1- 
parameter semigroup e*^. Although the continuous case is technically more difficult 
than the discrete cases, the results are extremely similar. Indeed, if one has a 
sequence of rational numbers Pi/qi converging to the irrational number z, then 
the statistical properties of Til{z) may be obtained as limits of the corresponding 
properties of Til{pi/qi). (Note that the reverse does not hold. One cannot obtain 
the statistics of a rational tiling by taking a limit of irrational tilings.) 

As always, we consider a basic right triangle Tq with sides a and b and hy- 
potenuse c. Let a = ln(c/a), (3 = ln(2c/6), and assume that z = a/ (3 is irrational. 
By a triangle of size s, we mean a triangle, similar to Tq, with hypotenuse ce~^ . 
Note that larger values of s correspond to smaller triangles, just as in the ratio- 
nal case, where the size Dk of triangles decreased with k. In our tiling the size 
parameter s will take values in [0, /u), where fu. = max(a;,/3). 

We now describe a semigroup similar to $'^. Let 5" be a collection of tiles, 
all with size in [0,/u). Expand this collection by a linear factor e*, resulting in 
triangles with sizes in [— t, // — t). Then subdivide the largest triangle, subdivide 
the largest remaining triangle, and so on, until all triangles have non-negative size 
parameter. By Lemmas 1 and 2, this occurs in a finite number of steps, and results 
in a collection of tiles with sizes in [0, fx). This collection is e^^S. 
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The semigroup e*^ naturally acts on the distribution of sizes. The properties 
of this action are summarized by the following theorem, which should be compared 
with Theorem 6: 

Theorem 8 

1) The eigenvalues X of L are the roots of p{X) = 0, where 

p(A) = e^^ - e^^-")^ - 4e(^-'')^. (5.1) 

There are no multiple eigenvalues. 

2) X = 2 is an eigenvalue. All other eigenvalues have real part strictly less than 
2. If a < (3, all eigenvalues have real part greater or equal to the real root of 
qI3X _|_g(/3-«)A _4^ while if a > /3, all eigenvalues have real part greater or equal 
to the real root of e"^ + 4e("~^)'^ - 1. 

3) The eigenfunction i{j{s), for fixed eigenvalue X, take the form 

^(s) = e^' - e^("-")/i(s -a)- 4e^^'-^^h{s - /3), (5.2) 

where h{x) is the Heavy side function 

h{x) = l^ ^/a^>0; ^53^ 

1 otherwise. 

4) Given an interval I C [0, n), the number of tiles with size in I is asymptotically 
a fraction 

\ f ds{a^h{a - s) + b^h{(3 - s))e'^' (5.4) 

of the total. 

5) Given an interval I C [0,//), the area covered by tiles with size in I is asymp- 
totically a fraction 

-^ -— I dsa^hia- s)+b'^hil3- s) (5.5) 

a^a + b^p J I 

of the total. 

Corollary: Let T be a tile of size s. Let Nt{s) be the number of tiles in e^^{T). 

Then 

2c^e~'^'' 
lim e-2*iVi(s) = (5.6) 

t^oo a?a + b^(3 

Proof of Corollary: From the distribution (5.4), we compute the average area per 
tile to be ab{a^a + 6^/3) /4c^. A tile of size s, rescaled by a factor e*, has area 
^^g2t-2s^2. Dividing by the area per unit tile we obtain (5.6). | 
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Proof of Theorem 8: For < t < min(a, /?, |a — /3|), the action of e*^ on population 
distribution functions is 

(iP{s + t)+A^{s + t-(3) if sG [/3-t,/3); 
{e^^^){s) = l'il){s + t)+i){s + t-a) ifsG[a-t,a); (5.7) 

(^ V('S + ^) ciU other s E [0, /u), 

where i^is) and e^^'(p{s) are understood to be zero for s > /U or s < 0. e*^ acts 
continuously on functions that are continuous away from 0, a, (3, and for which the 
oscillations at a and (3 are given by 

^{a-) - V'(a+) = /(O); ^{r) - V'(/9+) = 4/(0). (5.8) 

Setting 6*^-0 = e*'^V ^^ see that the eigenfunction ^{s) must equal e'^* times a 
piecewise constant function with discontinuities at 0, a, /?. Applying the boundary 
conditions (5.8), we obtain the eigenfunction (5.2). For s > /U, i^{s) then equals 
g(s-M)-'^ times p{\). The vanishing of ^(s) for s > // is equivalent to the eigenvalue 
equation p{\) = 0. Thus eigenfunctions satisfying the boundary conditions are in 
1-1 correspondence with roots oi p{X), with the eigenfunctions given by (5.2). 

Suppose a < /3, in which case p takes the form 

p(A) =e^^-e(^-")^-4. (5.9) 

A = 2 is a root, since p{2) = e'^'^ - e^^^^"") -4 = ^-^-4 = 0. This is the 
only real root, insofar as p{X) is an increasing function of A for A > 0, and p(A) is 
negative for A < 0. In particular, p{X) < implies that A < 2. 

Now consider complex roots X = Xr + iXj. If A/ 7^ 0, e^^ and e*^^~"-*^ cannot 
both be real, insofar as (3 is not a rational multiple of (3 — a. Thus |e^'^|, |e'^~"-*^| 
and 4 satisfy a strict triangle inequality. In particular, 

> le^^l - |e(^-")^| - 4 = p{Xr), (5.10) 

so Afj < 2. Also, 

< le^^l + |e(^-")^| - 4 = e^^« + e^^-")^^ - 4, (5.11) 

so Ai^ is greater than the real root of e^^ + g(/3-«)-'^ _ 4_ 

Now we exclude the possibility of multiple roots. A multiple root would require 
p(X) = p'{X) = 0. Suppose = p'{X) = /3e^^ + (/3-a)e(^-")^. Then e^^ and e^^"")^ 
must have the same phase, and their difference must also have that phase. However, 
if A is not real, e^^ and e^^~'^''^ cannot both be real, so their difference is not real, 
so their difference is not 4. Thus p'{X) = implies that p{X) 7^ 0, and there are no 
multiple roots away from the real axis. On the real axis, the only root is A = 2, and 
we have already seen that p'{2) is positive, not zero. 
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This establishes statements 1-3 for the case a < /3. The argument for a > P 
is almost identical, and is not repeated. 

Now suppose that we initially have a population distribution that is a linear 
combination of the eigenfunctions ip\. Applying e*^ to the system, the total area 
grows as e^*, so the number of tiles is bounded, both above and below, by a multiple 
of e^*. Applying e*^ and dividing by the number of tiles damps out all the modes 
with eigenvalue less than 2, i.e. all eigenvectors other than t(j2- In this case, the 
final distribution of sizes approaches a multiple of 'i/'2 5 and statements 4 and 5 of 
the theorem follow. 

Unfortunately, we cannot a priori assume that the initial condition is a linear 
combination of eigenfunctions of L, or that a test function is a linear combination 
of eigenfunctions of the dual operator L*. e*^ and e*^ are neither finite-rank 
operators nor self-adjoint operators on a Hilbert space, so standard theorems about 
the completeness of a basis is eigenfunctions cannot be applied. In principle it is 
possible for a test function / to have the property that e^^^ "^V does not converge 
at all. We must show that, when f = Xi (the characteristic function of an interval 
/), e^^^ "^V c?oes converge. 

Lemma 9: Given an interval I G [0, /u). The fraction of the area o/e*^(To) covered 
by tiles with size in I approaches a limit as t ^ oo. 

Given this lemma, it follows that the distribution of area of 6*^(5*), for any 
collection of tiles S, approaches a limit, from which it follows that the distribution 
of population of e^^{S) also approaches a limit. Since these limits are invariant 
under further evolution, and since the total area is proportional to e^*, these limits 
must correspond to the A = 2 eigenvector of L, hence must take the form (5.4) and 
(5.5). I 

Proof of Lemma 9: Assume a < (3; the other case is similar. Let Fi{t) be the 
fraction of area of e*^(To) covered by tiles with size in /. Given an e > 0, we will 
show how to compute a number such that, for all t sufficiently large, Fi{t) is within 
e of this number. Since this can be done for any e, limt^oo Fi{t) must exist. 

It suffices to show that the eventual fraction in an interval of size A, entirely in 
[0, a) or in [a, /3), and with A sufficiently small, can be estimated to within O(A^). 
Any larger interval can be broken up into a finite number of such small pieces, such 
that Yl errors < e. So let us fix an interval /, centered at sq, with width A. 

The strategy is this: We begin with an exact expression for the fraction of 
area of e*^(To) represented by tiles of size s. We sum this over s G / to get an 
exact formula for Fj{t). By taking certain limits and replacing certain sums with 
integrals, we obtain an expression that is independent of t. In the process we 
introduce two types of errors. One type can be made arbitrarily small by requiring 
t to be sufficiently large. The other type is O(A^). 

How many triangles of size s appear in e^^{To)7 That depends on whether s + t 
can be written as nia + n2(3 for non-negative ni and n2- If s -|- 1 = nia + n2/3, then 
a triangle of size s may be obtained by taking a triangle of size —t, subdividing 
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it, picking a daughter, subdividing it, picking a daughter, and so on for ni + 712 
subdivisions, with the descent involving ni daughters of type A and n2 daughters of 
type B. If s > a, then the last daughter must be of type B, or else after ni +77-2 — 1 
steps we would have already obtained a tile of size in [0,/9), and would not have 
made the final subdivision. If s < a there is no such constraint. We thus have 

r ("i+'^2^4'^2 [fs + t = nia + n2l3 and s G [0, a); 

Number of tiles of size s = I (n^'+n.-i^n, [f g + t = ma + n2l3 and s e [a, (3)- 

K otherwise. 

(5.12) 

Since at each division a fraction a^ jc? of the area goes into the A daughter, while 
a fraction 6^/c^ goes into the B daughters, then the fraction of the total area 
represented by tiles of size s is 

("'nT) (f^j \-^) if s + t = nia + mH and s G [0, a); 

rrr') (S)" (S)" if « + t = ni« + n2/? and. G [«,/?); ^^'^^^ 
otherwise. 

Note that 77-2 = (s + t — q;?ii)//3. Now let 

((.....(.-.))/.) /^r /gy---")/^ if.G[0,a); 



(s+t+n(/3-a)-/3)//3Waf ^ / .,2 \ ^ ^ > 1 1^ 



and let Jp be the periodic J-function 

5p(x) = ^5(x-n). (5.15) 



ne^ 



Note that /(s, t, n) is well-defined even when (s + 1 + n{(5 — a))/ 13 is not an integer. 
For t large, /(s, t, n) is a slowly-varying function of s and n. 

We then compute, exactly, 

F,(t)= J]sG/F,(t)= 5^ f{s,t,n)Sp{{s + t-na)/l3) (5.16) 

n=0 •^■f 

Next we approximate, by replacing /(s, t, n) by /(sq, t, n), where sq is the midpoint 
of /. This introduces an error that is a fraction 0(A) of the total. Since the total 
will turn out to be 0(A), the error introduced is O(A^). We thus have 

n=0 "^-f 
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Next we use the fact that multiples of an irrational number are uniformly 
distributed on M/Z. For any desired degree of accuracy, one can find an A^ such 
that, for any N consecutive integers n^, fj Sp{{s + t — nio)/ (3) equals one a fraction 
A//3 of the time (to within the allowed error), and the rest of the time. If we 
restrict ourselves to t so large that f{so,t,n) is nearly constant as n varies over 
intervals of size N, then 

Y, /(«o, t, n) / 6piis + t - na)/(3) ^ (A//3) J^ /(«o, t, n). (5.18) 

n=0 -^^ n=0 

Finally, the t ^ oo limit of the resulting sum can be expressed as an integral, and 
yields a nonzero number. | 

We now turn, as in §4, to the joint distribution of sizes and orientation. We 
parametrize the space of possible orientations of a single size, as in §4, by two copies 
of the unit circle, with Haar measure dfl. (See Figure 8). The space of all possible 
tiles up to translation, which we denote X, is 2 copies of 5"^ x [0, fx). 

Let / be a function on X . Given a collection 5" of tiles, let < /, 5" > be the 
average value of / on the individual tiles in 5", where each tile is given equal weight. 
Let < f, S >' be the average value of / on the individual tiles in 5", where each tile 
is weighted proportionally to its area. Let (/, 5") be the sum of / on the individual 
tiles of S. We define a semigroup e*^ acting on functions on X by 

(e*^7,^) = (/,e*^(5)). (5.19) 

Theorem 9 Assume a tiling Til{z), with z irrational and with 9 = tan~-'^(a/6) 
an irrational multiple of it. Let f be a continuous function on X , and let {Sn} 
he a sequence of supertiles of increasing size. Let du = e^^{a^h{a — s) + b'^h{(3 — 
s))dVtds/2c^ and dp = 2c^ e~'^^ dp / [a^ a + 6^/3) he measures on X , where h is the 
Heavyside function (5.3). Then 

lim < /, Sn >= I fdv (5.20a) 



>oo 



X 



hm < /, Sn >'= / fdp. (5.206) 

These limits also apply if f is the characteristic function of a rectangle in X . 

Remark: The measures du and dp are closely related to the integrands in (5.4) 
and (5.5), respectively. Theorem 9 states that the joint distribution of sizes and 
orientations is a product: the size distribution previously found in Theorem 8 times 
a uniform distribution of orientations. 

Proof: The proof is extremely similar to the proof of Theorem 7. As in that case, 
it is sufficient to establish (5.20a) for an arbitrary continuous function /. Such a 
function can be written as an absolutely convergent sum of Fourier modes (with 
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respect to rotations). The coefficient of eacli mode is a C^ valued function of s. 
The operator e*^ commutes with rotation, and so acts separately on each Fourier 
mode. 

On the n-th Fourier mode, e*^ , for t small, acts as follows. 

tL*A(._SAij{s-t + a) + Bij{s-t + (3) ifsG[0,t); . . 

■'J^' l^(s-t) another sG [0,;u), ^ ' 

where the matrices A and B are, as in §4, 

/ expHn(^ + 7r/2)]\ 

^" Vexp[m(^ + 7r/2)] J' ^^"^^^ 

r^ _ /exp(m^) + exp(m(^ + 7r)) 2exp(-m^) \ (koq) 

\ 2ex.p{in6) ex.p{—in6)+exp{in{—6 + 7j:))J' 

The 0-th Fourier mode decouples into (-^) and (_-^) components. The (-^) 
component is the distribution of sizes regardless of orientation, and its asymptotic 
behavior was already computed in Theorem 8. We must show that the n = (j^) 
component, and all the Fourier modes with n 7^ 0, grow strictly slower than the size 
of the system, and so represent a decreasing fraction of the system. 

We will control the L^ norms of the unwanted Fourier modes. To do this we 
need the L^ norms of the matrices A and B, and various products of A and B. 
The L^ norm of a matrix is maximum, over all columns, of the sum of the absolute 
values of the entries in that column. One can get a bound on the growth of the L^ 
norm of a mode by the mode with its absolute value, and replacing the matrices A 
and B by their norms. 

For n = 0, -B(_]^) = 0. With S = it is as if there is only one daughter per 
division, hence the L^ norm of the (j^) mode at time t is bounded by the L^ norm 
of the mode at time 0. Hence, as a fraction of the system, this mode shrinks like 

Next we consider n odd, for which the diagonal terms in B vanish. The sum 
of the absolute values of the entries of each column of B equals two. This is as if, 
at each subdivision, only two daughter B tiles are produced, instead of 4. To put 
it another way, at each subdivision a fraction 6^/2c^ of the area is lost. Since each 
piece of a tile of size —t must be divided at least t/(3 times, this means that the L^ 
norm of the n-th mode, for n odd, can grow no faster than e^*[l — (6^/2c^)]*'^ and 
so, as a fraction of the system, goes to zero. 

Finally we consider n even but nonzero. Here the column sums of A and B are 
the same as in the n = case, namely 1 and 4, respectively. However, the L^ norms 
of various products of A and B are smaller that in the n = mode. For example, 
B^ = 4cos{n9)B has norm 16| cos(n^)|, which is strictly smaller than 16. The norm 
oi BAB is also 16| cos(n^)|. Indeed, the only words in A and B which have norms as 
large as in the n = case are A"^ and A^'^BA^'^. Since the expansion of e*^ /, for 
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t > P, involves expressions such as B^, the growth of the L^ norm of n-th Fourier 
mode is bounded by an exponent strictly less than 2. As a fraction of the system, 
the n-th mode goes to zero. | 

Remark: The spectrum of e*^ may be obtained exactly as in §4. In seeking 
eigenvectors, equation (4.20) is replaced by 



(e^^ - e^^-^^^S - e^^-")^A)i; = 0. 



(5.24) 



Subsequent analysis may be repeated word for word, replacing A by e^, p by a, q 

by /9, m by /U, ifjk by ij{s), and r'^ by e~*. 

§6. Two exceptional tilings — Ti/(l/3) and Til{2) 

We saw in §4 how the population statistics of a rational tiling Til{p/q) depends 
on p and q. If g > 1, the second eigenvalue of the population matrix is greater than 
one, and fluctuations increase with subdivision. This leads to phenomena such 
as slippage along fault lines and a failure to be globally edge-to-edge. A typical 
example, Ti/(l/2), was studied in §3. 

In this section we study examples of the remaining cases. We study Til{2) as 
an example of a Til{p/1) tiling. In all such tilings, the second eigenvalue of the 
population matrix is less than one. We shall see how, in the case of Ti/(2), this 
prevents slippage along fault lines. 




Figure 10a 




Figure 10b 




Figure 10c 



Figure 10. Three stages of subdivision for Til{2). 

Finally, we consider Ti/(l/3), the only rational tiling to exhibit only a flnite 
number of orientations of each size of tile. As in all cases with q > 1, there is 
an eigenvalue greater than one in the problem. The fluctuations governed by this 
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eigenvalue are enough to force Tz/(l/3) to be nonperiodic. However, because of 
rational relations between the lengths of certain edges, the tiles in Til{l/3) meet in 
only a finite number of ways. 

We begin with Til{2). Til{2) is based on the right triangle with a = v5 — 2 ^ 
0.2361, b — 2vv^ — 2 ^ 0.9717, c = 1. Several iterations of the subdivision are 
shown in Figure 10. An essential feature of Til (2) is 

Theorem 10 The triangles in the tiling Til{2) meet in only a finite number of 
ways. 

Proof: The proof is essentially in two steps. First we show that the slippage along 
the primary fault line is bounded. Then we show that slippage along a fault line is 
the only means by which tiles in a rational tiling can meet in an infinite number of 
ways, and that all fault lines are similar to the primary fault line. 

To examine what happens along the primary fault line, we consider the bound- 
ary of T2n, the 2n-th subdivision of the basic triangle. Note that the hypotenuse 
and short leg of T2 consist only of hypotenuses and short legs of big triangles. Ap- 
plying the subdivision again, we get that the hypotenuse and short leg of T4 also 
consists only of hypotenuses and short legs of big triangles. Similarly for all T2n- 

As in the proof of Theorem 5, we consider the evolution of these legs and 
hypotenuses as a one dimensional substitution system in its own right. Let a denote 
the effect of subdividing the basic triangle twice. Under a, each hypotenuse (denoted 
H) is replaced by 4 hypotenuses and a short leg (denoted 5"), while each short leg 
is replaced by a hypotenuse. That is, the one dimensional population matrix is 

M=(^ I), (6.1) 

with eigenvalues A± = 2 ± y^ and eigenvectors v± = { J . 

Now let Hn and 5"^ be the number of hypotenuses and short legs in a^{H), and 
let H'^ and S'^ be the number of hypotenuses and short legs in a'^{S). By expanding 
(-1^) and (q) in terms of v±, it is easy to see that 

{V5-2)H^-S^ = -{2-^/5r+'; {Vh - 2)^- S'^ = -{2 - VhT . (6.2) 

Next we measure slippage. Let P and R be the endpoints of a hypotenuse, as 
in Figure 7, and let E be any intermediate point, not necessarily the midpoint. Let 
fn{E) be the number of complete short legs, between P and E, in a^{PR), minus 
the number of short legs between P and E in a^{RP). As in the proof of Theorem 
5, fn{E) measures the extent to which the two tiles of T2n+2 that meet at E are 
offset. 

Lemma 10 \fn{E)\ < 5. 

Proof: Let vq = P, and let Vk be the vertex of a^{PR), between P and E, that 
is closest to E. Note that, in a^[PR), there are at most 4 hypotenuses and at 

31 



most one short leg between Vk-i and Vk, since the interval Vk-iVk was only part 
of a hypotenuse or short leg in a'^~^{PR). By (6.2), (y^ — 2) times the number of 
hypotenuses in a'^{PR) between Vk-i and Vk, minus the number of short legs, is 
bounded in absolute value by {\/5 — 2)"^"'^. Summing over /c, we get that {\^ — 2) 
times the number of complete hypotenuses between P and E in a^{PR), minus the 
number of short legs, is bounded in absolute value by X]i^o(^~^)* < 2. A similar 
bound applies to the number of hypotenuses and legs in a'^{RP). Thus the surplus 
of short legs on one side of PE relative to the other, plus (\/5 — 2) times the deficit 
of hypotenuses, is bounded by 2 + 2 = 4. Since a surplus of short legs implies a 
deficit of hypotenuses, the surplus of short legs is itself bounded by 4. | 

We return to the proof of Theorem 10. Lemma 10 limits the number of ways 
for two tiles to meet across a fault line. Suppose two tiles ti and t2 meet across 
a fault line PR, modeled on a'^{H). Let E he a, point on their common edge. 
The distance from the vertex of ti closest to P to the vertex of t2 closest to P is 
either \fn{E)\{^/5 — 2) or 1 — \fn{E)\{^/5 — 2). Since |/„| is at most 4, this means 
there are only a finite number of ways for two triangles to meet across such fault 
lines. Now a^{H) is the result of subdividing the basic triangle an even number of 
times. However, since subdivision is deterministic, having only a finite number of 
distinct configurations in the even subdivisions implies that there are only a finite 
number of distinct configurations in the odd subdivisions, and thus a finite number 
of configurations in all. 

To complete the proof of Theorem 10, we must show that every pair of adjacent 
triangles either meets full-face to full-face, or meets across a fault line based on 
successive subdivision of a hypotenuse. Consider two tiles, ti and ^2, that meet. 
Let k be the smallest integer such that both tiles lie in the same supertile S of order 
k. Since ti and t2 do not meet in a supertile of order k — 1, ti and t2 must meet 
across one of the five lines of the first subdivision of 5". See Figure 10a. There is 
a local refiection symmetry across edges 1 and 3, so if ti and t2 meet across these 
edges they must meed full-face to full-face. Edge 2 is the primary fault line. Further 
division (see Figure 10c) shows that there is local refiection symmetry across edge 
4, while edge 5 is a hypotenuse-based fault line, as considered above. | 

We now turn to the tiling Tz/(l/3). By Theorem 4, Ti/(l/3) is the only tiling 
in our construction to have both a finite number of sizes of tiles, each of which 
appears only in a finite number of orientations. Ti/(l/3) is based on an isosceles 
right triangle. Although the two legs have the same length, we distinguish between 
the two, calling the "6" side "long" and the "a" side "short", in analogy to the 
tilings with b > a. In subdividing we must specify which legs of the daughter tiles 
are labeled "long" and "short". This is shown in Figure 11, and several further 
subdivisions are shown in Figure 12. 
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Figure 11. "Long" and "short" edges in Ti/(l/3). 






Figure 12. Three stages of subdivision for Tz/(l/3). 

Theorem 11 The tiling Tz/(l/3) is nonperiodic. The tiles meet in only a finite 
number of ways. 

Proof: As usual, we consider the 1-dimensional substitution scheme induced on the 
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edges by 2-fold substitution. The 2n-th subdivision of the hypotenuse consists of 
some hypotenuses of large triangles (denoted if), some long legs of medium sized 
triangles (L), and some hypotenuses of small triangles (h). Let a denote the action 
of two subdivisions, a takes each H to an H and two L's, each L to two /I's, and 
each h to an. H. The population matrix is 



M= 2 (6.3) 




with eigenvalues 2 and (— 1± V^7)/2, and with the eigenvector | 1 | corresponding 

to eigenvalue 2. Note that the complex eigenvalues have magnitude v2. 

Asymptotically, the three types of edges appear in a ratio of 1:1:1. However, 

since the initial condition I is not an eigenvector, and since all eigenvalues are 

greater than 1 (in magnitude), the difference in number between if s and L's, or 
L's and /I's, with grow exponentially with time. That is, while the total population 
grows as 2"^, the fluctuations grow as ^/2 . 

If Tz/(l/3) were periodic, a long line of the form a^{H) would consist of several 
periods, plus a remainder at each end. Each period would have H^s, L's, and /I's in 
exactly a 1:1:1 ratio, so only the partial periods at each end could contribute to the 
difference in population between H and L. Thus the population difference would 
remain bounded as n ^ oo. Since this difference is unbounded, Ti/(l/3) cannot be 
periodic. 

Finally we note that the length of h is the same as that of L, and half that 
of H. This simple ratio of lengths means that slippage along the fault line has no 
effect on the number of ways triangles can meet. If two tiles meet across a" (if), 
either they have a vertex in common or their closest vertices are separated by the 
length of h. Thus there are only a finite number of ways for two tiles to meet across 
a fault line. By the same argument as in the proof of Theorem 10, this implies that 
there are only a finite number of ways for triangles to meet at all. | 

§7. Conclusions 

We have constructed a family of substitution systems, indexed by the parameter 
z = ln(sin(^))/ ln[cos(^)/2], where 9 is an angle in the basic triangular tile. We have 
established the following properties. 

1) The tilings generated by these substitutions are all non-periodic. 

2) The tilings have well-defined limiting distributions of size and orientation. If 
O/tt is irrational, this distribution is rotationally invariant. In Radin's termi- 
nology, the tilings have "statistical rotational symmetry". The form of the 

34 



joint distribution of size and orientation suggests that the tihng has a purely 
absolutely continuous spectrum. 

3) The tilings with rational z all satisfy the hypotheses of Goodman- Strauss 's 
theorem, implying that they can be forced through local matching rules. 

4) The rational tilings Til{p/q), with q > 1, have statistical fluctuations that 
grow with iterations of the substitution rule (although they grow slower than 
the size of the system). In Til{l/2) these fluctuations force triangles to meet in 
an infinite number of distinct ways. We conjecture that this infinite diversity 
of local behavior is a property of all rational tilings Til{p/q) with q > 1 and 
p/q ^ 1/3. 

5) In the rational tilings Til{p/1) the eigenvalues that control fiuctuations are all 
less than one. In Til{2) this forces the tiles to meet in only a finite number of 
local patterns. We conjecture that this is a property of all tilings Til{p/1). 
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